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We show that a combination of linear absorption spectroscopy, hyper-Rayleigh scattering, and a 
theoretical analysis using sum rules to reduce the size of the parameter space leads to a prediction 
of the two-photon absorption cross-section of the dye AF455 that agrees with two-photon absorp- 
tion spectroscopy. Our procedure, which demands self-consistency between several measurement 
techniques and does not use adjustable parameters, provides a means for determining transition 
moments between the dominant excited states based strictly on experimental characterization. This 
is made possible by our new approach that uses sum rules and molecular symmetry to rigorously 
reduce the number of required physical quantities. 

PACS numbers: 42.65.An, 33.15.Kr, 11.55.Hx, 32.70.Cs 



I. INTRODUCTION 

There is a long history of using nonlinear-optical tech- 
niques to build an understanding of the mechanisms 
of light-matter interactions. Given the availability of 
mostly single- wavelength lasers, early measurements used 
time domain studies to deconvolute mechanisms such as 
molecular reorientation and the electronic response in 
liquidsQI, 0, H 3 which were used to make fast optical 
gates, [5| and in solids, for example, to study excitons in 
quasi-one-dimensional polymeric crystals. @, 0] 

One of the first attempts to use dispersion in the Op- 
tical Kerr Effect (OKE) to understand the nonlinear- 
optical response in liquids was based on a qualita- 
tive comparison of experimental results with the sum- 
over-states (SOS) expression for the nonlinear-optical 
susceptibilities [8j as calculated by Orr and Ward.[9[ More 
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recently, OKE dispersion measurements and more sophis- 
ticated multi-state models have been applied to deter- 
mining the transition moments between excited states in 
silicon phthalocyanine-monomethacrylate. 

The difficulty with such approaches is that they re- 
quire either non-realistically simple models with only a 
couple of parameters; or, complex models that involve 
transition moments between excited states that can not 
be independently verified by experiment. Furthermore, a 
set of parameters that successfully models one particular 
measurement is not often consistent with other indepen- 
dent measurements. 

The sum rules are quantum mechanical identities that 
relate the dipole matrix elements and energies to each 
other; so, the SOS hyperpolarizability can be ex pre ssed 
in terms of a subset of the dipole matrix. [Ill, [T3, 13] In- 
deed, the dipole-free SOS expression is derived by using 
the sum rules to eliminate all terms with dipole moment 
differences. Dipole-free expressions for both the hyper- 
polarizability, /3,[H| and the second hyperpolarizability, 
7, [la] have been derived and shown to be mathematically 
equivalent to the standard SOS results. These dipole-free 
expressions are key to significantly reducing the number 
of parameters required to model the nonlinear response. 

In the present work, we use linear absorption spec- 
troscopy to determine the energies and transition mo- 
ments from the ground state of the octupolar chro- 
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mophore commonly referred to as AF-455, and use hyper- 
Rayleigh scattering, the symmetry of the molecule, 
and the dipole-free SOS expression for the hyperpolar- 
izability to determine the transition moment between 
the two dominant states. Using no adjustable parame- 
ters, we predict a two-photon absorption spectrum that 
agrees with experimental results. This suggests that our 
approach may be a simple alternative that is straight- 
forward to apply and yields self-consistent results that 
spans linear absorption, hyper-Rayleigh scattering, and 
two-photon absorption spectroscopy. 

What makes our work unique is that all of the quan- 
tities required to predict the linear and nonlinear-optical 
response are determined experimentally, which is made 
possible using sum rules to reduce the number of re- 
quired parameters. Other approaches have been intro- 
duced that provide estimates of the nonlinear response 
from simple measurements. For example, Perez-Moreno 
and coworkers have introduced a rule of thumb that pro- 
vides a rough estimate of the resonant two-photon ab- 
sorption (TPA) cross-section of a large set of molecules 
simply using the number of it electrons. [16] Rebane and 
coworkers used a theoretical approach based on the den- 
sity matrix to show that parameters determined from 
a measurement of the linear absorption spectrum of a 
dipole transition can be used to determine the TPA cross- 
section at the one-photon absorption maximum. [17] For a 
broad range of molecules, the approach was shown to be 
yield TPA cross-sections that deviated at most by 50% 
compared with measurements. However, while these ap- 
proximate techniques are useful for estimating the TPA 
cross-section, they do not predict the dispersion of the 
TPA cross-section, nor do they predict the first hyper- 
polarizability. In contrast, the approach presented here 
leads to an accurate prediction of the TPA spectrum. 



II. THEORY 

Our approach begins by simplifying the analysis of the 
dispersion of the nonlinear-optical susceptibilities by pro- 
viding a model that depends only on a reduced set of mea- 
surable molecular parameters. In contrast, most studies 
reported in the literature rely on calculating and/or us- 
ing as fit parameters transition energies between excited 
states, excited state dipole moments, and energies. 

Because excited-state parameters can not be experi- 
mentally verified, and often many sets of parameters can 
yield reasonably good fits of the data to the theory, the 
conclusions based on such studies may not be sound. Fur- 
thermore, one set of parameters will often provide a good 
fit to one experiment, but not to others. To compound 
such problems, semi-empirical calculations of a nonlinear 
susceptibility measurement off-resonance does not cor- 
rectly predict the resonant behavior. 

We begin by introducing the SOS dipole-free expres- 
sions of the nonlinear-optical susceptibilities. These ex- 
pressions reduce the number of parameters needed to de- 



scribe the nonlinear-optical response by eliminating the 
dipole terms. Then, we use the symmetries of an oc- 
tupolar molecule to show that the two-fold degenerate 
first excited state can be expressed in a basis where the 
transition moment between these two states vanishes. 
However, this basis demands that the dipole moments 
be non-vanishing. This is of no consequence since we 
use the dipole-free expressions to calculate the nonlinear 
susceptibilities. Finally, applying the sum rules to the 
dipole basis, we show that the four- level system can be 
expressed as a three-state model. 



A. Sum Rules and the Dipole-Free SOS 
Expressions 

The polarizability of a molecule along x for an incident 
field of frequency u> polarized along x is given by,@ 
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where x nm is the transition moment between states n and 
ra,[32| E h q — E n — Eq is the energy difference between 
states, and where the prime denotes the fact that the 
ground state is excluded from the sum. The energies are 
complex to account for damping. Note that while the 
susceptibilities are tensors, for the sake of simplicity, we 
express only the largest tensor component. 

The sum rules are calculated from the Schrodinger 
Equation using the closure identity. For charges under 
the influence of electric and magnetic fields, the sum rules 
are exactly obeyed, and are derived from the matrix ele- 
ments of commutators of the Hamiltonian H and position 
x, which yields [l~2| 
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Note that the set of relationships given by Equation [2] 
with I — p are commonly called the sum rules. The more 
general form above with I ^ p was developed to account 
for the off-diagonal components. [ill [l2| Equation [5] cor- 
responds to an infinite set of equations, which are each 
labeled by distinct pairs of integers I, p. The sum rules 
clearly show how the matrix elements of the position op- 
erator (proportional to the dipole matrix) and energy 
levels are intimately related to each other and therefore 
can not be independently adjusted. 

The infinite set of equations embedded in the sum 
rules can be used to simplify the SOS expressions for 
the nonlinear-optical susceptibilities. The off-diagonal 
sum rules can be used to express the ground and excited 
state dipole moments in terms of transition moments. [l4| 
As such, the dipole terms can be eliminated. [3] Using 
these relationships, the common SOS expression can be 
transformed into the dipole-free SOS expression, which 
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FIG. 1: The three degenerate states of an octupole with no 
charge hopping between sites. 



for the largest diagonal component of the hyperpolariz- 
ability yields, [3 
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and where u> m rj = w° n0 — «7 m o. fiuJ,Q is the energy differ- 
ence between excited state |m) and the ground state and 
7 m o is half the natural linewidth for a transition from 
state m to the ground state. This expression is some- 
times called the reduced SOS expression. 

The second hyperpolarizability, as well as higher-order 
hyperpolarizabilities, can be transformed into a dipole- 
free form. The resulting algebraic expressions are too 
complex to present here, but their form can be found in 
the literature, flal 



B. The Dipole Basis 

We begin by using the symmetries of an octupolar 
molecule with D^h symmetry to define the energy level 
diagram. It is typical to model an octupole with a two- 
fold degenerate excited state. This degeneracy can be 
understood by considering a basis in which the charge is 
placed on one of the three branches of the octupole as 
shown schematically in Figure [1] With no charge hop- 
ping, the states |1), |2), and |3) are the three eigenstates 
of the system. 

Using the approach described by Fevnman.fl9| we de- 
fine a 120° rotation by the operator R. Because a rotation 



of the system by 360° leaves the system invariant, R 3 = 1, 
and the eigenvalues of R are clearly 1, exp(27ri/3), and 
exp(47r«/3). With hopping between nearest neighbors, 
the Hamiltonian is, 



H = 1£ -AR- AR- 
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where £q is the energy without hopping and A the hop- 
ping energy. Since H and R commute, they share the 
same eigenfunctions, which are of the form: 

M = ^(|l) + |2) + |3)), 

h/>i> = ^=(|l) + e 2 "/ 3 |2) +e 4 "/ 3 |3)), 



and 



h/< 2 } = i=(|l)+ e 4 W3| 2)+e 2W3| 3) ). (6) 



y/3 

The corresponding Eigenenergies are, 

E = £ -2A 
and 

E x = E 2 =£ + A. 



(7) 



These three states describe the ground state and the 
two-fold degenerate excited state of the octupole. If the 
excitation on branch n designated by the state vector |n) 
represents a dipole of magnitude |/x| = e- 1, where I is the 
distance between the center of the molecule and the end 
of one of the branches, then the dipole moment of state 
i is given by, 



fi = (tpi\ r\i>i 
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where the unit vectors 1, 2, and 3 are along the three 
branches of the octupole represented by the state vectors 
|1), |2), and |3), respectively, as diagrammed in Figure 
[TJ and, — e • I ■ n is the dipole moment of the branch 
represented by the state vector \n). Note that we have 
used 
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under the assumption that there is no overlap between 
the wavefunctions on different branches. Because each 
state has no dipole moment, we call this the octupolar 
basis. 

Since the three-level model may not be sufficient to 
describe an octupole, we add a third excited state. The 
properties of this state are completely general and there 
are no restrictions on its dipole moment. Figure [2] shows 
the energy level diagram of the four-state system that we 
will use for our analysis and labels all of the elements of 
the dipole matrix. Note that the degenerate states are re- 
labeled 1 and 1' to stress that these states are degenerate 
with energy E\ . 
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FIG. 2: The four-level model of an octupolar molecule. All 
possible transition moments along the i-direction are shown. 



In the octupolar basis, r*oo: rn, and fyy vanish. It is 
simple to show that r±y ^ 0, and will be of the form, 
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Since any linear combination of two degenerate energy 
eigenvectors is an energy eigenvector, it is possible to 
chose linear combinations of |</>i) and 1^) to form two 
orthogonal energy eigenvector with x ■ r\y — 0; but then, 
rn and ryy will no longer vanish. In this new basis, the 
vectors in the degenerate subspace are given by, 

I'M = ^(l^» + (l^» = ^(2|l>-|2)-|3)). 

I0i'> = - (IV>2» = ^= (|2) - |3)) . (11) 

We call this the dipole basis. Note that \<f>o) — |?/>o) and 

|0 2 ) = |^a). 

The dipole moments of state \4>i) and \4>y) in the dipole 
basis are then given by, 

1 
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where we have used the coordinate convention shown in 
Figure[T] Clearly, the dipoles in the two degenerate states 
point in opposite directions. It is straightforward to show 
that the transition moment r\y has no i-component in 
the dipole basis: 



1 



ni' oc 



-2 + 3) oc y. 



(13) 



Finally, in the dipole basis the transition dipole mo- 
ment rot' has not i-component: 
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(14) 



C. Applying the sum rules to the dipole basis 

Our approach is to determine the hyperpolarizability 
and second hyperpolarizability along x, that is, to cal- 
culate (3 XXX and j xxxx , and then to use octupolar sym- 
metry to determine the other components. As such, we 
need only consider the one-dimensional sum rules. Fur- 
thermore, we will assume that four states are sufficient to 
describe the nonlinear response, and will therefore only 
consider four-state sum rules. 

We evaluate the sum rules in the dipole basis, which 
is chosen deliberately to meet the condition x\y = 0. 
Equation 1141 yields x$y = 0. To summarize, the dipole 
basis demands: (1) the no coupling condition, x\y = 0; 
(2) the transition moment from the ground to primed 
state vanishes, Xoy — 0; and (3) the ground state mo- 
ment vanishes, xoo = 0. We note that the dipole basis 
does not impose any other restrictions on the moments. 
However, as we show below, the sum rules impose an 
additional set of conditions. 

The sum rule I = and p = V yields, 

- X 00 X yE W + X Q1 XiyE W + X yXyyE W 

+ XQ2X2V {2E 20 - E w ) = 0, (15) 

which in the dipole basis (with xq§ — 0, x§y — and 
x\y = 0) yields 



xqiXiv (2£ 20 _ -Em) = 0. 
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Since (2E 2 q — Eiq) > 0, the implication is that X02 or 
X2i> must vanish. We choose x 2 i> = because it is in 
agreement with the sum rule I = 1 and p = 1', which 
yields 

- X 10 X yE 1Q + XuXxyEu 

+ X 1 yXyyE 11 +X 12 X 2 yE 21 = 0. (17) 

As a consequence of the fact that En = 0, 

£12*21' (£20 - £10) = 0. (18) 

Since E20 — E\q > 0, we conclude that x 2 y = 0, which is 
in agreement with Equation 1161 

The net result is that all transitions that include state 
1' are forbidden. We call this channel blocking. As a 
result of channel blocking, the four-level system can be 
described by a three-level model. Note that we could 
have alternatively made the choice that state 1 is the 
excluded state (i.e. 1-^1'). This relabeling is equivalent 
to rotating the coordinate system by 90°. 

We note that while the choice of one basis over an- 
other results in different dipole moments and energies of 
the degenerate states, each basis will yield the same ob- 
servables - such as the absorption spectrum, hyperpolar- 
izability, etc. So, we choose the dipole basis because it is 
the most convenient: the number of parameters required 
to describe the nonlinear response is reduced and the 
dipole-free form of the nonlinear susceptibilities makes 
the values of the ground and excited state dipole mo- 
ments irrelevant. 



D. The Polarizability 

In the four-state model, the sum rule with I = and 
p = yields, 
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Since the transition moment to the second excited state 
and the energies of the degenerate states and second ex- 
cited state are independent of the basis, 
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The polarizability near resonance with the degenerate 
state (E\o ~ fuv) is give by Equation [T] 
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The area under the peak in the absorption spectrum 
near E\q = huj is proportional to |a:oi| 2 + |xoi' | 2 , and no 
measurement can separate the individual contributions 
of each degenerate state. Since the sum rules demand 
that any choice of basis in the degenerate subspace must 
yield |xoi| 2 + | a^oi' 1 2 = constant, all bases will lead to 
the same observed linear absorption spectrum. Thus, 
the choice of basis has no effect on any observable. The 
dipole basis is the particular choice where xov — so all 
of the oscillator strength in the subspace is concentrated 
in scoij and, the area under the linear absorption peak 
centered at E\q will be proportional to | £Cqi | 2 Ew- 



E. The Hyperpolarizability in the Dipole Basis 

In this section, we impose octupolar D^h symmetry 
and the sum rules to calculate the hyperpolarizability 
using the dipole-free expression in the dipolar basis. In 
analogy to the rigorous proof that the polarizability is 
independent of basis, the same will hold for the hyperpo- 
larizability, which we state here without proof. 

Recall that the dipole-free SOS expression for the hy- 
perpolarizability has numerators of the form xo n x nm x m o 
where m ^ n and m/0 for all m and n. It is convenient 
to represent all possible virtual transitions graphically, as 
shown in Figure [3] 

In the dipole basis, diagrams (e) and (f) are forbidden 
because of the no coupling condition. Diagrams (c) and 
(d) are forbidden by channel blocking 

The only non- vanishing terms represented by diagrams 
(a) and (b) represent the three-level dipole-free SOS ex- 
pression for [3 XXX . We note that the numerators associ- 
ated with these diagrams are complex conjugates of each 
other. The contribution to the hyperpolarizability of one 
diagram is the same as the other one with states 1 and 2 
interchanged. 
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FIG. 3: All possible contributions to the dipole-free form of 
the hyperpolarizability in a four-state system. 



The Second Hyperpolarizability in the Dipole 
Basis 



In this section, we impose octupolar D$h symmetry 
and the sum rules to calculate the second hyperpolar- 
izability using the dipole-free expression in the dipo- 
lar basis. The dipole-free SOS expression for the sec- 
ond hyperpolarizability has two types of terms. In the 
first type, called Type I, the numerators are of the form 
XQnXnmXmixio where m ^ n and m ^ I. Also, no index 
can represent the ground state, 0. It is convenient to 
represent all possible virtual transitions graphically, as 
shown in Figure 2] 

Diagrams (c) through (f) are forbidden by channel 
blocking. Diagrams (g) through (1) are forbidden be- 
cause of the no-coupling condition. As in the case of the 
hyperpolarizability, the non-vanishing diagrams (a) and 
(b) represent the two terms in the three-level dipole-free 
SOS expression for the second hyperpolarizability. 

Next we consider Type II terms, as shown in Figure [5] 
Diagrams (b) and (c) are forbidden by channel blocking. 
The only non- vanishing term is once again the three-level 
dipole-free Type II term. 



III. RESULTS AND DISCUSSION 

The dipole basis, when used with the dipole-free SOS 
expression, leads to a model of the linear and nonlin- 
ear susceptibilities that depends only on the energies of 
the two excited states E\q, E20; their widths; and the 
transition moments £01, ^02, and x\2- The quantities 
E10, E20, xqi, and 202 can be determined using the lin- 
ear absorption spectrum. The only remaining quantity 
is x\2- The transition moment between the two excited 
states can be determined from a measurement of /3 at 
one wavelength provided that the other quantities are 
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FIG. 4: All possible contributions to the dipole-free form of 
the second hyperpolarizability in a four-state system of the 
Type I terms. 
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FIG. 5: All possible contributions to the dipole-free form of 
the second hyperpolarizability in a four-state system of the 
Type II terms. 




FIG. 6: The AF455 molecule. Note its D 3h symmetry. 
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known a priori. From these measured quantities, we will 
show that all-optical and nonlinear-optical properties are 
predicted with no adjustable parameters. 

We apply this technique to the octupolar molecule AF- 
455 in solution, which as shown in Figure [6j has D^h 
symmetry. The absorption spectrum is used to deter- 
mine its polarizability using the fact that with appropri- 
ate choice of molecular coordinate system, a xx — a yy 
and a xy = a yx = 0. Our analysis holds generally for 
any octupolar molecule with D^h symmetry, which we 
will refer to simply as an octupole. An isotropic solu- 
tion of octupoles will have an isotropic polarizability of 
(a) = 2a xx /3. Thus, from an absorption spectrum, the 
polarizability along x can be directly determined. 

The linear absorption spectrum of a sample is obtained 
by measuring the transmittance of a broad spectrum 



FIG. 7: Extinction spectrum of AF455 and fit to four Gaus- 
sians; three to fit the first excited state and one for the second 
excited state. 



source passing through a cuvette of 1cm path length, 
and containing a solution of the sample in solvent. The 
spectrum is referenced to the transmittance through pure 
solvent in an identical cuvette. Figure [7] shows the mea- 
sured extinction spectrum of AF455 in tetrahydrofuran 
(THF). 

The two peak positions are used to determine the en- 
ergies E\a and E20 ; and , the width at half maximum for 
each peak is used to determine r 10 and r 2 r). The absorp- 
tion spectrum is fit to four Gaussian peaks; three for the 
first excited state and one for the second excited state 
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as shown in Figure [7l The sum of the areas of the first 
three peaks is used to calculate the transition moment 
to the first excited state, using the method described by 
Tripathi et al,[2(| while the area of the fourth Gaussian is 
used to get the transition moment to the second excited 
state. The solid curve is the sum of all four peaks, and 
shows that the four-Gaussian theory provides a good fit 
to the data. The results are summarized in Table |TJ 

The hyperpolarizability is determined through hyper- 
Ray leigh scattering (HRS), using a femtosecond pulsed 
laser (Millennia X + Tsunami with a lock-to-clock sys- 
tem that ensures 80Mhz pulsed output) at fundamental 
wavelength of A = 800nm in conjugation with a low fre- 
quency lock in amplifier and a signal generator. Details 
can be found elsewhere. [Hi [HJ The measured HRS hy- 
perpolarizability yields, 
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(22) 



The octupolar symmetry of the compound (which is used 
in Equation 1221 to determine (3 XXX directly from an HRS 
measurement) was confirmed by the depolarization mea- 
surements. Note that the other components of the (3 
tensor are given by,[23[ 



ftxxy 



ftxy 



ftyxx i 
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while all other tensor components vanish. A demodu- 
lation technique is used to determine a fluorescence-free 
value to insure that only the hyperpolarizability is be- 
ing measured. [24] We note that for AF455, no fluores- 
cence contribution was found at 400nm, which leads to 
a more accurate determination of the first hyperpolariz- 
ability than is possible with the demodulation technique 
when fluorescence is present. 

The dipole-free expression for the rra;iE-component of 
the first hyperpolarizability for a three-state model is 
given bv[14j 
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where D~^(o;i, UI2) is a dispersion function that depends 
only on the energies of the first two excited states, the 
two photon frequencies, and the widths of the two states. 
In the HRS measurement, u)\ — 0J2 = uj. We note that 
for an infinite number of states, the standard SOS and 
dipole-free expressions are rigorously identical; but when 
truncated to three states, the two are different. (Note 
that our sum rule/symmetry analysis represents a four- 
level system with a three-level model without loss of in- 
formation.) It is not possible a priori to know which 
expression is more accurate. Our strategy is to use the 
dipole-free expression because it does not require knowl- 
edge of the dipole moments. The merits of this approach 
will be judged by the predictive capability of the model. 
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FIG. 8: The experiment used for measuring the two-photon 
absorption cross-section. 



The linear absorption measurement determines all the 
transition moments, energies and widths except for Xyx- 
Thus, given that the first hyperpolarizability is measured 
at a known frequency ui, Equation 1241 can be inverted to 
solve for /ii 2 = — ex ±2- Table U shows all of the parame- 
ters determined from linear absorption spectroscopy and 
one measurement of the first hyperpolarizability. Note 
that the transition dipole moment is related to the posi- 
tion matrix elements, jj, nm = — ex nm . The three values 
of /ii2 listed in Table U represent the uncertainty range 
due to uncertainty in the HRS measurement. 

Figure [5] shows the experiment used for determining 
the two-photon absorption cross-section using the mea- 
sured two-photon fluorescence power. This technique was 
developed by Xu and Web. (2f| Details of how the data is 
related to the two photon cross-section can be found in 
the original paper. [25| The advantage of this experiment 
is that it is a reliable method for determining the TPA 
cross-section[26| and is not as susceptible to excited state 
absorption as is nonlinear transmission. [2?], [28| Here we 
briefly describe those issues that are particular to our 
implementation of the technique. 

The sample solution is prepared by adding 0.0958 gram 
of AF455 crystals as received from Wright Patterson Air 
Force Base to 200ml optical spectrum grade THF, in a 



Fitting Method E 10 eV E 20 (eV) | Ml0 |(D) l/j 20 |(D) r 10 (eV) r 20 (eV) |Mi2| Upper (D) |^i2| mlddle (D) |mi2 | lower (D) 
3 Gaussian 3.0 4.1 12.6 9.4 0.22 0.35 15.5 12.6 9.7 



TABLE I: Experimental parameters (from the linear absorption spectrum) and |/ii2| measured using HRS. The range of HRS 
values reflect experimental uncertainties. All quantities are dressed values. 



clean flask at room temperature. The mixture is agi- 
tated in an ultrasonic water bath for thirty minutes to 
make a uniform solution. Two quartz cuvettes (ordered 
as a matched pair) are filled, and labeled SI and S2, re- 
spectively, are filled to 4/5 full with the uniform solution. 
These two identical samples are used to calibrate the col- 
lection efficiency of the reference and sample arms of the 
TPA experiment. 

The reference solution is made from 200ml of 100M 
Rhodamine B solution by adding 0.00958 grams of Rho- 
daminc B powder to 200ml optical spectrum grade 
methanol followed by agitation in an ultrasound bath for 
30 minutes at room temperature. A clean cuvette is filled 
with this solution to 4/5 full and capped tightly. 

A sample is excited with a laser beam, and the two- 
photon fluorescence signal (integrated over all fluores- 
cence wavelengths) is detected with an RCA C31034A- 
02 photomultiplier tube (PMT). This PMT is of high 
quantum efficiency over a broad range of wavelengths 
throughout the visible. Filters are used to remove any 
wavelengths corresponding to the pump energy or less, 
leaving only the fluorescence spectrum due to two-photon 
absorption. The time-integrated current from the PMT 
is proportional to the number of two-photon absorptions. 
A typical run uses the average over many lasers pulses of 
the 10 Hz laser source to decrease noise. The reference is 
used to take into account laser fluctuations and is used 
as a standard for determining the absolute two-photon 
absorption cross-section. 

For molecules of D^h octupolar symmetry, the tensor 
components are related to each other according to.[23j 

Jxxxx — 'fyyyy — ^Jxxyy — ^Jxyxy — ^tC . . . (^^) 

Thus, knowledge of the value of j X xxx allows one to de- 
termine all other components using Equation 1251 We use 
the dipole-free expression to calculate Jxxxx because it 
does not require knowledge of the ground and excited 
state dipole moments. [15| We do not present the dipole- 
free-expression here, but rather refer the reader to the 
literature. fl5| 

Figure [5] shows the theoretically-calculated value of the 
imaginary part of j X xxx &s a function of wavelength us- 
ing the measured values of the dipole moment matrix 
elements, energies, and widths shown in Table |TJ The 
upper and lower curves show the uncertainty range due 
to the experiential uncertainty in determining the excited 
state transition moment fj,\2 from HRS. 

The two-photon absorption spectrum, as measured 
with the two-photon fluorescence experiment, yields an 
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FIG. 9: The measured two-photon absorption spectrum 
(points) and the theoretical curve (middle), where all quan- 
tities are converted into the imaginary part of jxxxx- The 
upper and lower curves represent the uncertainty introduced 
due to the uncertainty in the HRS-determined transition mo- 
ment /ii2. Note that there are no adjustable parameters in 
the theory. 

isotropic average over all tensor components of the two- 
photon absorption spectrum. Given the D^h symmetry 
of the molecule, the isotropic value can be related to the 
individual tensor components. Thus, the measured TPA 
cross-section can be related directly to the imaginary part 
of Jxxxx ■ The measured values are shown as points with 
error bars in Figure [U We note that since the samples 
used in all measurements are liquid solutions, all quanti- 
ties such as the polarizability, hyperpolarizability, transi- 
tion dipole moments, etc. are dressed values. [29| Vacuum 
quantities can be determined using the appropriate local 
field models. H| 

The theoretical spectrum and the data are in good 
agreement and all but one of the data points fall within 
the error band of the predicted spectrum. Thus, the ap- 
proach of using the dipole-free expressions for the first 
and second hyperpolarizability in the dipole basis pro- 
vides a theoretical description that is consistent with 
three sets of measurements. In particular, the predicted 
TPA spectrum through the two-photon resonance is con- 
sistent with the data. This suggests that our approach 
may be useful in modeling the dispersion of the linear 
and nonlinear susceptibilities of octupolar molecules with 
D 3h symmetry. 
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In light of the fact AF455 is a complex molecule, it 
may appear somewhat surprising that the predicted dis- 
persion of the TPA cross-section using a four-level model 
is in such good agreement with the data. This may be due 
to two factors. First, the HRS measurement was deter- 
mined near two-photon resonance, where the TPA peak 
is measured so the contributions of the first excited states 
are heavily weighted and dominate the response. In con- 
trast, an off-resonant HRS measurement potentially in- 
cludes contributions from the tails of many higher-energy 
excited states, thus yielding an inaccurate determination 
of the transition moment /ii2- On-resonance measure- 
ment insures that the influence of the transition moment 
fii2 is large. Secondly, since the TPA spectrum is mea- 
sured only near the two-photon resonance, the same set 
of states are being probed. Finally, it is possible that 
the observed agreement is a coincidence. Similar studies 
of other octupolar systems would determine the general 
applicability of our method. 

It is instructive to apply the same approach to the 
three-level model using the standard SOS expressions 
under the assumption that the three states are non- 
degenerate and have no dipole moment. Using the value 
of Xi2, determined form HRS measurements and the 
standard SOS expression form the hyperpolarizability, 
the theoretical value of the second hyperpolarizability 
predicts a two-photon absorption spectrum that is two- 
orders of magnitude larger than the measured one. This 
illustrates how the typical approach can lead to inaccu- 
rate results, and supports the validity of our approach. 

IV. CONCLUSION 

We have introduced an approach that combines mea- 
surements, symmetries, and sum rules to fully character- 
ize the important states of a molecule that allows all of 
the linear and nonlinear susceptibilities to be accurately 
modeled with no adjustable parameters. Our approach 
is general in that it can be applied to any molecule of 
any symmetry class. In the present work, we have illus- 
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